Magnetic long-range order induced by quantum relaxation in single-molecule magnets 
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Can magnetic interactions between single-molecule magnets (SMMs) in a crystal establish long- 
range magnetic order at low temperatures deep in the quantum regime, where the only electron 
spin-fluctuations are due to incoherent magnetic quantum tunneling (MQT)? Put inversely: can 
MQT provide the temperature dependent fluctuations needed to destroy the ordered state above 
some finite Tc, although it should basically itself be a T- independent process? Our experiments on 
two novel Mn4 SMMs provide a positive answer to the above, showing at the same time that MQT 
in the SMMs has to involve spin-lattice coupling at a relaxation rate equaling that predicted and 
observed recently for nuclear spin-mediated quantum relaxation. 

PACS numbers: 75.40.Cx, 75.45.-|-j, 75.50.Xx 



Despite the large number of studies on magnetic quan- 
tum tunneling (MQT) in molecular crystals of single- 
molecule magnets (SMMs) (ll, the question whether it 
is able to bring the spin system into thermal equilib- 
rium with the lattice, remains unsolved. Prokof'ev and 
Stamp have suggested that interaction with rapidly 
fluctuating hyperfine fields can bring a significant num- 
ber of electron spins into resonance. Coupling to a nu- 
clear spin bath indeed allows ground-state tunneling over 
a range of local bias fields ^ much larger than the tunnel 
splitting A, whereas, in its absence, tunneling would hap- 
pen only if ^ is A or less. Within this theory, magnetic 
relaxation could thus in principle occur with no exchange 
of energy with the phonons of the molecular crystal [Sj . 
Support for the Prokof 'ev/Stamp model came from mag- 
netic relaxation studies on the Feg SMM However, 
these experiments covered only initial stages of the re- 
laxation process, leaving open the question whether the 
final state corresponds to a thermal equilibrium. When 
MQT is combined with coupling to a heat bath, dipo- 
lar couplings between cluster spins can induce long-range 
magnetic ordering (LRMO) at sufficiently low tempera- 
tures 0. This phenomenon has not been observed yet 
for any of the known SMMs relaxing by MQT. Attempts 
have been made in Feg and Mni2 in particular, by increas- 
ing the tunneling rate by means of an applied transverse 
field Q. However, the failure is inherent to the approach, 
since the applied fields needed are much higher than the 
weak interaction energies involved, and will thus destroy 
the LRMO. 

In this letter, we present zero-field time-dependent spe- 
cific heat measurements performed on two novel tetranu- 
clear molecular clusters, both with net cluster spin S — 
9/2, denoted by Mn4Cl and Mn4Me, which have similar 
cluster cores but different ligand molecules. For both, the 
cluster spins become blocked along their anisotropy axes 
for temperatures in the liquid helium range, and relax- 



ation below ~ 0.8 K can thus only proceed by incoherent 
MQT between the two lowest lying states m = ±5'. For 
both compounds we prove below that the MQT has to 
be inelastic. For Mn4Me, the tunneling rates are found 
sufficiently high to establish, even at zero field, ther- 
mal equilibrium conditions down to the lowest tempera- 
tures (~ 0.1 K). Accordingly, the MQT channel enables 
the occurrence of LRMO between the cluster spins at 
Tc — 0.21 K. Comparing the magnitude of Tc with Monte 
Carlo simulations suggests the coexistence of dipolar and 
weak superexchange interactions between clusters. In 
view of the essential role of the dynamic nuclear bias in 
the MQT mechanism, our results call for an extension of 
the nuclear-spin mediated quantum relaxation model 01 
to include inelastic processes, where MQT is accompa- 
nied by phonon creation or annihilation. 

Analytically pure samples of the compounds 
Mn403L(dbm)3, with L = Cl(0Ac)3 or [02C(C6H4- 
p-Me)]4, hereafter abbreviated as Mn4Cl and Mn4Me, 
were prepared as described in Ref. 0- AH samples were 
characterized by elemental analysis. Both molecules 
possess a distorted cubane core with one Mn''^ ion 
(spin S — 3/2) and three Mn'^+ ions (S = 2), superex- 
change coupled via three oxygen ions. The mtra-cluster 
exchange couplings have been studied by magnetic 
susceptibility measurements Below T < 10 K, the 
four Mn spins become ordered with a net total spin 
S = 9/2 subject to an uniaxial crystal field 8], the 
symmetry axis running approximately through the Mn"*"*" 
ion and the L- ligand. Whereas Mn4Cl has a local virtual 
Csv symmetry, the more bulky carboxylate ligand 
distorts and lowers the symmetry (Cg) of Mn4Me in 
the crystal and increases the magnitude of the uniaxial 
anisotropy 

Low-temperature specific heat measurements were per- 
formed by a home-made calorimeter using a thermal re- 
laxation method, see Ref. P|. By varying the thermal 
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resistance of the thermal Hnk between calorimeter and 
cold-sink, the characteristic timescale Te of the experi- 
ment can be varied. In this way time-dependent specific 
heat measurements can be exploited to investigate spin- 
lattice relaxation when the relaxation time becomes of 
the order of Te (~ 1 — 100 s) 0]. The samples consisted of 
1 — 3 mg of polycrystalline material, mixed with 2 — 5 mg 
of Apiezon-N grease to ensure good thermal contact. The 
high-temperature (2 K < T < 300 K) specific heat was 
measured for a Mn4Cl pellet sample of about 40 mg using 
a commercial (PPMS) calorimeter. 

The specific heat C/R of Mn4Cl is shown in Fig. 1, 
which combines data obtained for Te « 2 s, 8 s, and 300 s 
(as estimated at T = 0.4 K) with the ones obtained with 
the high-T calorimeter. Let us first consider data mea- 
sured for T > 1 K, where C is independent of Tg. A 
A- type anomaly is observed at T ~ 7 K, having a relative 
height of 2.5 R. Measurements with an applied magnetic 
field (not presented here) have shown that this anomaly 
is insensitive to the field, proving that it has to be as- 
sociated to a structural phase transition. Between 1 K 
and 7 K, C is dominated by contributions from the lat- 
tice phonons and from transitions between energy levels 
of the 5 = 9/2 multiplet split by the crystal field. 

Neglecting, in first approximation, the Zeeman split- 
tings ^ due to effective fields representing the hyperfine 
interaction and the mter-cluster magnetic couplings, the 
spectrum of energy levels is doubly degenerate at zero- 
applied-field, so that C only depends on transitions be- 
tween levels inside each of the potential wells that are 
separated by the anisotropy energy barrier U . We shall 
call this the intra-weW contribution. The associated mul- 
tilevel Schottky anomaly is calculated with the eigenval- 
ues of the spin Hamiltonian 

n = -DSl - E{Sl - SD + AiSt (1) 

with parameters D, E and obtained independently 
from inelastic neutron scattering and high-frequency 
EPR measurements ^| . The intra- well specific heat 
decreases exponentially as T decreases and, since the first 
excited m = ±7/2 level is about {2S — 1)D = 5.5 K above 
the TO = ±9/2 ground-state doublet, it becomes almost 
neghgible when T < 0.8 K (Fig. 1). Adding the lat- 
tice contribution, calculated with a Debye temperature 
~ 15 K, to the Schottky accounts well for the exper- 
imental data (solid line in Fig. 1). The lattice specific 
heat above 10 K appears to be composed of a number of 
Einstein- type contributions (not shown). 

Below 1 K, we expect the equilibrium magnetic spe- 
cific heat {Cm) to be dominated by two contributions. 
The first arises from incoherent MQT events inside the 
ground-state doublet that is split by the action of the 
effective fields arising from hyperfine interactions and 
mter-cluster dipolar coupling. The second is the specific 
heat Cnuci of the nuclear spins of Mn, whose energy lev- 
els are split by the hyperfine interaction with the atomic 
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FIG. 1: Temperature dependence of the specific heat of 
Mn4Cl measured for Te ~ 2 s, 8 s and 300 s. Solid line, 
sum of the Debye term of the lattice (dotted line) plus the 
Schottky contributions; dashed line, calculated nuclear con- 
tribution C„ucl- 

electron spins. The dashed line in Fig. 1 represents Cnuci 
calculated with the hyperfine constants A^f = 7.6 mK 
and Ahf — 11.4 mK for, resp., Mn'^+ and Mn^+ ions, 
obtained from ESR 12]. Experiments performed for the 
longest Te ~ 300 s show indeed a large low-T contribu- 
tion. By contrast, the specific heat decreases by almost 
two orders of magnitude when Te « 2 s, evidencing that 
Te has a large effect in this temperature range. This 
shows that the equilibrium between the relative popula- 
tions of the ±9/2 and —9/2 states cannot be established 
within Te if this is too short. We note that, for the short- 
est Te, the low-T specific heat becomes even smaller than 
Cnuci , indicating that both nuclear and electron spins are 
out of equilibrium. This is understandable, since the only 
channel for the nuclear spins to exchange energy with the 
lattice is via the electron spins. The strong connection 
between nuclear and electron spin-lattice relaxation has 
earlier been observed for Mni2, Mug and Feg [^Il3|. 

The experimental C/R of Mn4Me, measured for Te ~ 
4 s (as estimated at T = 0.4 K), is shown in Fig. 2 
together with the calculated contributions of the lat- 
tice and the nuclear spins, and the intra-weW Schot- 
tky contribution. The experimental data display a A- 
anomaly at Tc ~ 0.21 K that we attribute to the on- 
set of LRMO. The Schottky anomaly is calculated as 
for MmCl with parameters obtained from high-frequency 
EPR liallll. For T > 1 K, the remaining contribution 
is given by the lattice and is well described by the sum 
of a Debye term for the acoustic low-energy modes plus 
an Einstein term for a higher energy mode, with val- 
ues 9d — 12.3 K and 6*^; ~ 22 K for, resp., the Debye 
and Einstein temperatures. Below 0.15 K, the specific 
heat of Mn4Me shows a clear upturn that can be de- 
scribed by Cnuci /R ^ 4.27 x lO'^/T^. This contribution 
is well fitted using hyperfine constants Ahf = 8.7 mK and 
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FIG. 2: T-dependent specific heat of Mn4Me measured for 
Te « 4 s. Daslied line, calculated nuclear contribution Cnuci', 
solid line, sum of lattice (dotted line) plus Schottky contribu- 
tions. Insets: electronic Cm and entropy variation AS, after 
subtraction of Cnuci- Solid line, calculated dipolar ordering 
for the easy axis along the (110) direction (see text); dashed 
lines, high-T entropy limits for 5* — 1/2 and S = 9/2. 

Ahf = 13.8 mK for, resp., Mn^+ and Mn''+ ions. These 
values are of the same order as for Mn4Cl and are close 
to those reported by Zengh and Dismukes for a natural 
Mn4 cluster 

As seen in Fig. 2, the magnetic ordering peak and the 
calculated Schottky due to the splitting of the S — 9/2 
multiplet are well separated in temperature. This is 
confirmed by the analysis of the temperature depen- 
dence of the electronic entropy AS{T) = J^[Cm{T) — 
Cnuci{T)]/TdT (right inset of Fig. 2). As expected, the 
total entropy of the electron spins tends to i?ln(2S' -t- 1), 
with S = 9/2, at high temperatures. However, for the 
magnetic ordering region {T < 0.8 K), AS corresponds 
to an effective spin S = 1/2, as appropriate for a two- 
level system. This proves that there only the two lowest 
levels (m — ±9/2) are populated and contribute to Cm- 

We note that, although the anisotropy barrier U ~ 
DS^ - AiS^ -\E\S^ amounts to w 14 K for both Mn4Me 
and Mn4Cl, due to the lower symmetry of Mn4Me the 
2nd-order off-diagonal term (E) in the Hamiltonian of 
Eq. (0), and thus also A of the ground-state, will be 
much larger for this compound. This is confirmed by 
high-frequency EPR experiments that give EJD ~ 0.21, 
i.e., nearly 5 times larger than for Mn4Cl pH llll|. We 
observe that this difference has a very large influence on 
the spin-lattice relaxation. For similar Te values, the elec- 
tron spins of the Mn4Cl molecule go off equilibrium below 
T = 0.8 K (Fig. 1), whereas for Mn4Me we observe ther- 
mal equilibrium for electron and nuclear spins down to 
the lowest temperature. We conclude this from the fact 
that (j) the total electronic entropy contribution equals 
the expected limit for S = 9/2; (ii) the remaining specific 
heat below 0.15 K agrees well with the calculated Cnuci- 



As commonly found in molecular clusters, the mag- 
netic core of Mn4Me is surrounded by a shell of non- 
magnetic ligand molecules. It follows that, mter-cluster 
superexchange interactions are very weak, leaving the 
mier-cluster dipolar coupling responsible for magnetic 
ordering To check this for Mn4Me, we have per- 

formed Monte Carlo simulations, as described in Ref. j^, 
for a S — 9/2 Ising model of magnetic dipoles regu- 
larly arranged on the Mn4Me lattice. We have repeated 
our calculations for several orientations of the molecular 
easy axis (an example is given in the left inset of Fig. 2 
for the easy axis along the (110) direction). The calcu- 
lated Tc's are always much smaller than the experimen- 
tal value. Consequently, the Mn4Me molecules are also 
coupled by weak superexchange interactions. Indeed, by 
adding an mier-cluster nearest-neighbor exchange inter- 
action \J\/kB — 0.14 K to our dipolar calculations, we 
reproduce the experimental Tc value Q. 

In order to estimate the spin-lattice relaxation rate 
r for Mn4Cl at low temperatures, we have used the 
relation for the time-dependent specific heat Cm{i) ~ 
Co+{Ceq — Co) [l — e~'"'^'=] , where Co and Ceq are the adi- 
abatic and equilibrium limits of the specific heat, respec- 
tively 0]. For the electron spins, Co is to good approx- 
imation given by the m^ra- well Schottky specific heat, 
whereas the "slow" specific heat at equilibrium corre- 
sponds to excitations involving transitions between the 
two wells. We have fitted, thus, the Cm(Te) data of 
Mn4Cl, taking at each temperature for Ceq the Cm{T) 
measured for Mn4Me (in the range T > Tc)- We show 
r(T) in Fig. 3, together with data from ac-susceptib ility 
and (short-time) magnetic relaxation experiments |l5| . 
For T > 1.7 K, r follows the Arrhenius law, with 
U ~ 13.5 K and To « 1.4 x 10"^ s. For T < 0.8 K, 
r deviates from this thermal-activation law, in remark- 
able agreement with the magnetic relaxation data. Its 
weak temperature dependence confirms that relaxation 
to thermal equilibrium is dominated by direct MQT tran- 
sitions within the ground-state doublet. 

Summing up, we have proven that for both compounds 
thermal contact between spins and lattice is established 
by MQT fluctuations. For Mn4Me, the associated rate is 
even fast enough to produce thermal equilibrium down 
to the lowest temperature and thus enable LRMO. This 
implies F > I/te ~ 1 s~^. For Mn4Cl, we find a lower 
rate {10~^ - 10"^ s''^) for T < 0.8 K. It is of interest to 
compare the experimental F for Mn4Cl with predictions 
from conventional models for spin-lattice relaxation, as- 
suming that the m = ±9/2 energy levels of the cluster 
spins are time-independent. We have simulated the ef- 
fect of mter-cluster dipolar coupling and hyperfine inter- 
actions by introducing static magnetic fields = 150 G 
and Bz — 350 G [l£j. The presence of these Zeeman 
terms is essential, otherwise tunneling would be forbid- 
den for half-integer spin S = 9/2 [13]. We calculated F 
by solving a master equation, including intra- as well as 
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FIG. 3: Spin-lattice relaxation rate of Mn4Cl: (*) and (•) [l^ 

obtained from ac-susceptibility and magnetic relaxation data; 
(o) obtained by fitting the Cm{Te) data of Fig. 1 (see text). 
Dashed line is fit of high-T data to Arrhenius law; solid line 
is calculated for magnetic fields Bx = 150 G and Bz = 350 G 
(see text). 



inter-well transitions, induced by phonons only, between 
exact eigenstates of the spin Hamiltonian of Eq. as 
discussed in Ref. and recently applied to the analysis 
of Cm of Fes and Mni2 clusters p. The result (solid 
line in Fig. 3) agrees well with the activated behavior 
observed at high-T, but fails to account for F measured 
below 1 K by six orders of magnitude! This large dis- 
crepancy can not be ascribed to errors in the estimated 
elastic properties of the lattice. Both the pre-factor of 
the Arrhenius law and the measured value of 0d give a 
value of Cs — 5(1) x 10^ m/s for the speed of sound. By 
contrast, F observed at low-T would require Cg and On to 
be 15 times smaller. This would give rise to a large lattice 
specific heat well below 1 K, which is not observed. 

It appears therefore that extension to dynamic hyper- 
fine fields acting on the cluster-spin levels, as proposed in 
Ref. |3| , is indeed a necessary prerequisite for any model 
for the MQT of SMMs. Such dynamic bias fields will 
sweep the tunneling levels with respect to one another, 
thereby enabling incoherent Landau-Zener type tunnel- 
ing events. The model predicts quantum relaxation rates 
agreeing with experiments 0, but so far the relaxation 
of the cluster spins was thought to occur solely /primarily 
to the nuclear spin bath, relaxation to phonons was only 
expected at much longer time-scales. Our specific heat 
experiments, in which obviously the heat is transferred 
to the spins via the lattice, clearly demonstrate that in 
fact spin-lattice relaxation has to be involved and at much 
the same fast rates! Since application of 'conventional' 
models for spin-lattice relaxation leads to rates orders of 
magnitude too low, our data call for an extension of the 
Prokof'ev/Stamp model in which nuclear spin-mediated 



MQT events are combined with creation or annihilation 
of phonons. 
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Stamp for enlightening discussions, and R.S. Edwards 
and S. Hill for HFEPR experiments on Mn4Me. This 
work is part of the research program of the Stichting 
voor Fundamenteel Onderzoek der Materie (FOM). 
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